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We present a double distribution function of dark matter halos, with respect to both object mass 
and local over- (or under-) density. This analytical tool provides a statistical treatment of the 
properties of matter surrounding collapsed objects, and can be used to study environmental effects 
on hierarchical structure formation. The size of the "local environment" of a collapsed object is 
denned to depend on the mass of the object. The Press-Schechter mass function is recovered by 
integration of our double distribution over the density contrast. We also present a detailed treatment 
of the evolution of overdensities and underdensities in fi m + Qa ~ 1 and Q m = 1 universes according 
to the spherical evolution model. We explicitly distinguish between true and linearly extrapolated 
overdensities and provide conversion relations between the two quantities. 
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I. INTRODUCTION 

Cosmological distributions have long been used with 
great success as an analytical tool complementary to nu- 
merical simulations. They have been used to constrain 
the cosmological parameters; interpret results of cosmo- 
logical simulations; study regions of the parameter space 
which cannot be approached by simulations due to pro- 
hibitive computational cost; exploring the effects of vari- 
ous physical processes in an efficient if approximate way. 
The analytical tool used most widely in cosmology is the 
mass function of dark matter halos (distribution of the 
number density of halos with respect to halo mass). Ana- 
lytical descriptions of dark matter halos are usually based 
on the Press-Schechter formalism ([1],[2]) and its exten- 
sions (e.g., [3], [4], [5], [6]). The Press-Schechter mass 
function has is in good agreement with the results of N- 
body simulations (e.g. [7], [8], [9]). More sophisticated 
approaches taking into account deviations from spherical 
symmetry (e.g., [10], [11], [12]) have improved this agree- 
ment even further. 

To derive the Press-Schechter mass function, regions 
in space are smoothed on successively smaller scales. 
The mass of a collapsed object is then taken to be the 
largest smoothing mass scale for which the average linear 
overdensity exceeds some threshold. In this way, mat- 
ter in the universe is distributed among collapsed struc- 
tures of different masses, which all share the same value 
of average overdensity (the threshold value). Informa- 
tion about the local environment of the collapsed objects 
(whether they live in underdensities or overdensities) is 
thus erased. 

For this reason, and despite its wide applicability, the 
expression for the mass function cannot be used to ad- 
dress environment-related questions: Does the mass func- 
tion of structures in superclusters differ from the mass 
function inside voids? Are structures of some particular 
mass more likely to reside inside underdense or overdense 
regions in space? How does such a preference evolve with 



redshift, and how sensitively does it depend on the cos- 
mological parameters? How does the state of the material 
surrounding and accreted by a collapsed object depend 
on the mass of the object and the cosmic epoch? 

To address such questions, we seek a double distribu- 
tion of the number density of structures with respect to 
mass but also to local overdensity (or underdensity) . In 
order to extract information about the surroundings of 
collapsed structures, we use the same random walk for- 
malism which rigorously yields the Press-Schechter mass 
function ([2], [3]). Integration of this distribution over 
density contrast should return the Press-Schechter mass 
function so that the successes of the Press-Schechter for- 
malism be retained. 

Two complications arise in the effort to expand the 
Press-Schechter mass function to incorporate a descrip- 
tion of the local overdensity. First, the concept of the 
"local environment" is somewhat vague and needs to be 
defined in a more rigorous way. The size of the "local en- 
vironment" cannot be the same for all structures. If this 
was the case, very small structures would represent only a 
tiny fraction of the "environment" , while very large struc- 
tures could even exceed the size of the "environment" , 
which would be an unphysical situation. This problem 
is not exclusive to analytical tools, but also needs to be 
addressed when analyzing the results of numerical sim- 
ulations. Second, the Press-Schechter treatment of the 
density field uses linear theory, and ways of converting 
this information to a more physical non-linear result need 
to be determined. 

We address the first problem by introducing a cluster- 
ing scale parameter, (3, which allows us to define the size 
of the "environment" of each structure as a function of its 
mass. We address the second concern by calculating con- 
version relations between the linear-theory overdensities 
(or underdensities) and those predicted by the spherical- 
evolution model. 

Environment-related questions in cosmological struc- 
ture formation have also beed addressed using analytical 
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models for the clustering properties of dark halos which 
evaluate quantities such as the cross-correlation function 
between dark halos and matter, and the biasing factor 
(e.g. [13], [14], [15], [16], [17], [18], [19], also see review 
by [20] and references therein) . These analyses are based 
on the same random walk formalism which we use here 
to derive our double distribution (also see [21] for fit- 
ting formulae from N-body simulation results). However, 
the information content of the double distribution, which 
treats the "environment" in a mass-dependent fashion, is 
complementary to that of correlation functions, which 
describe the clustering properties of the dark halo pop- 
ulation at some fixed spatial scale. The double distribu- 
tion is ideally fitted for population studies of cosmolog- 
ical objects. If the properties of a single object can be 
parametrized as a function of its mass and its environ- 
ment, then the double distribution can be used to predict 
the statistical properties of such objects, as well as their 
evolution with time, for any cosmological model. 

Our paper is organized as follows. In §11 we briefly re- 
view the "random walk" formalism used in the derivation 
of the Press-Schechter mass function by [2] and [3] . Wc 
then use the same formalism to derive the double distri- 
bution for dark matter halos. We also discuss converting 
between linear and spherical-evolution density contrasts, 
and we present some interesting derivative quantities of 
the double distribution. In §111 we explore the informa- 
tion content of our double-distribution by plotting the 
distribution itself as well as its derivative quantities for 
concordance (r2 m + CIa = 1) and Einstein-deSittcr uni- 
verses. Finally, we conclude and discuss our findings in 
§IV. 



II. FORMALISM 

A. Labeling of Cosmic Epochs 

Throughout this paper we use the value of the dimen- 
sionless scale factor of the universe, a, to label different 
cosmic epochs. We normalize it so that the scale fac- 
tor of the present is ao = 1. Then, a given value of a 
corresponds to a redshift 

z = a~ 1 -l, (1) 

independently of the cosmological model used. On the 
other hand, the conversion between a and time t does 
depend on the assumed cosmology. For a flat (fi to t = 
Ptot/pc = 1) universe containing only matter and vac- 
uum energy in the form of a cosmological constant (with 
present-day density parameters £! m = p m ,o/Pc,o an d 
= PA,o/y°c,o respectively), the conversion equation is 



(2) becomes 



t = -H^ni 1/2 sinh" 1 (V^aM 



(2) 



where Hq is the present-day value of the Hubble parame- 
ter. In the limiting case when = 0, O m = 1, equation 



t 



(3) 



We choose to cast our results in terms of a because it 
has a trivial and cosmology-independent relation with the 
directly observable redshift. In addition, it increases with 
time, which allows for a more intuitive interpretation of 
the evolution of the quantities we consider here. 



B. Random Walks and the Press-Schechter Mass 
Function 

The "random walk" formalism was introduced for the 
derivation of cosmological mass functions by [2] and by 
[3]. Here, we briefly review the basic concepts of this 
formalism, before extending it to the case of the double 
distribution in the next section. 

The Press-Schechter mass function of collapsed struc- 
tures is the comoving number density of virialized objects 
per differential mass interval, dn/dm, for every cosmic 
epoch a [32]. A related quantity is the mass fraction, 
P(> to, a), which is the fraction of matter in the uni- 
verse belonging to collapsed structures with mass > to. 
If P(> m, a) is known, then dn/dm can be calculated 
from 



dn p m0 
— (to, a) = — - 
dm to 



- — P(> to, a) 
dm 



(4) 



where p m .o is the present-day matter density of the uni- 
verse. 

P(> to, a) is in turn calculated by assigning, at each 
epoch a, every infinitesimal element dm in the universe 
to a collapsed structure of some mass to. A structure is 
considered "collapsed" if its mean overdensity 



(*> = 



(A 



'structure 



> - Pn 



(5) 



exceeds a certain critical value, S c (a). In equation (5), 
Pm,a is the mean matter density of the universe at epoch 
a. The critical oversendity S c (a) is the mean overdensity 
predicted by the spherical evolution model for a struc- 
ture virializing at epoch a. For each point in space, 
one calculates the mean local overdensity by smoothing 
the overdensity field S(x,a) with a spherically symmet- 
ric filter function of varying mass scale, starting from 
to — > oo , where one averages over the whole universe and 
finds identically (S) =0, and proceeding to successively 
smaller scales. When a mass scale is found for which 
the mean overdensity becomes equal to 6 c (a), it is taken 
to be the mass to of the parent object of the infinites- 
imal mass at the point under consideration. This way 
of assigning object masses circumvents the structure-in- 
structure problem, since the mass of the parent object is 
always the largest possible mass satisfying the criterion 
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for collapse. All information on substructure within col- 
lapsed structures is thus erased from the resulting mass 
function. 

The way the average overdensity (S) changes when 
the smoothing mass scale is varied resembles, under cer- 
tain conditions, a ID random walk [2]. For all "parti- 
cles" (points in space in our case), the walk begins at the 
"spatial origin" ((<5) = 0), at "time zero" (m — > oo). 
As "time progresses" (m decreases), each "particle" may 
move either to the "left" (negative (5)) or to the "right" 
(positive (5)). An "absorbing wall" exists at 8 c (a). If 
this "wall" is reached, the "particle" is "removed" from 
the walk (the point is assigned its parent object mass and 
removed from further consideration at smaller values of 
to). P(> m,a) is then the fraction of "particles" which 
have been "lost" by "time" m, and it can be calculated 
using random- walk theory. 

However, we must first ensure that simple random- 
walk theory is indeed applicable. First, each "step" of 
the "walk" should be completely independent from the 
previous step. This requires that the fc-modes producing 
an increase A (5) in the space-like variable not appear in 
any of the previous steps in (6). A smoothing window 
function sharp in k— space, 



W m (k) 



1 k < k c {m) 
k > kjm) 



(6) 



(see [2] and [3] for more extended discussions on the con- 
sequences of such a choice) enforces this condition, since 



and 



/ 



W„ 



\x -x\)S(x)d 3 x = (5 k e iSo i d 3 k (7) 

J k<k c (ro) 



A<5>* = / 

J k. 



5 k e lS °- k d 3 k 

k c (m)<k<k c (m — dm) 



(8) 



which only involves new fc-modes corresponding to scales 
from m to to — dm. 

Second, there must be an equal probability for the sys- 
tem to "move" towards any one of the two available "di- 
rections" . A Gaussian overdensity field (which is the 
usual assumption for deriving analytic mass functions 
and which we adopt in this paper) guarantees that this 
condition is satisfied. 

Finally, the appropriate "time-like" variable (which 
should depend on m) needs to be selected, given that 
the "space- like" variable is (S) . By direct analogy to the 
ID random walk theory result (a; 2 ) = 2Dt, and from the 
definition of the variance of the overdensity field S(m), 

r-k(m) 

S(m) = a 2 {m) = (\5{m, x)\ 2 ) = / k 2 dk\5 k \ 2 (9) 

Jk=o 

we can immediately identify Dt — > S(m)/2. 

Three further complications need to be addressed. 
First, our knowledge of 8 k and subsequently S(m) is 



limited at late times. In the early universe, right after 
matter-radiation equality, (|<5/c| 2 ) can be simply described 
in terms of a power-law in k modified by a transfer func- 
tion, (|4| 2 ) oc T 2 {k)k n . While all 6 are still in their 
linear regime, they simply grow by the linear growth fac- 
tor (independent of k). However, at later times, when 
certain structures start departing from the linear regime, 
we cannot use our simple early-universe expressions for 
5k- Second, (6) is limited to be > — 1, which introduces a 
second, reflecting "wall" at a value of (8) = — 1, further 
complicating the random- walk calculations. Finally, the 
true overdensity field loses its Gaussianity as it evolves 
past the linear regime. 

To circumvent these problems, we define the linearly 
extrapolated overdensity field, S(x,a), as the overdensity 
field that would result if all structures continued to grow 
according to the linear theory until time a. Now S(x, a) 
is not limited to be > —1, since it does not represent 
real overdensities. In addition, we can always calculate 
S(m) for S(x, a), since 5 k is modified from its simple 
early-universe expression only by the linear growth fac- 
tor. Finally, the extrapolated field remains Gaussian at 
all times. 

The linearly extrapolated overdensity 5(x, a) and the 
associated variance S(m), are time- varying, but the time 
dependence is well-known (see appendices B and C), and 
the same for both S and 8 2 [33]. Thus, the time depen- 
dence drops out of ratios 5/VS which appear in the mass 
function. For this reason we may, without loss of general- 
ity, choose any single epoch to evaluate these quantities, 
with the stipulation that 5 and S(m) must refer to the 
same epoch. Given that ct(to) is often normalized to the 
present value of erg, a convenient choice of epoch is the 
present. Then, equation (9) gives for S(m) 



S(m) 



, J k l7T 2 (k)k^ 2 dk 
' Xfio^ T 2 (k)k n + 2 dk 



(10) 



Thus we only consider 5(x, do) (the overdensity field 
linearly extrapolated to the present epoch), which we use 
instead of the true field 5(x, a) in our random walk for- 
malism [34]. To find the mass function at a particular 
cosmic epoch a, we calculate the location of the "absorb- 
ing wall", S c (a). If a structure is predicted to collapse 
at epoch a according to the spherical evolution model, 
then So, c (a) is the overdensity this same structure would 
have had if, instead of turning around and collapsing, it 
had continued its linear evolution until the present. This 
So,c{a) is then our "absorbing wall". 

We can now derive the mass fraction and mass function 
using random walk theory. If a particle executes a one- 
dimensional random walk with an absorbing boundary at 
a point xi, then its probability W(x,t) to be between x 
and x + dx at time t is [22] 



W(x, t, x\)dx = 



exp 



4Dt 



exp 



(2xi 



4Dt 



2V^Dt 



dx , 



(11) 
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where x < x\. In our case, the probability that a point 
in space will be assigned an average extrapolated over- 
density between 5 and S + dS when filtered at a scale to 
corresponding to a variance of S(m) is 



W(5,S,5 , c )d5 = 



exp 



il 

2S 



exp 



(2<5 , c -5) 2 



2S 



dS. 



(12) 

with S < (5o. c - The mass fraction P(> to, a) is then the 
fraction of points already "lost" from the walk when fil- 
tering at higher mass scales, which is one minus the frac- 
tion of points remaining in the walk, 



P(> to, a) = P(> 5 ,c) = 1-1 

J — ( 

= crfc ( ~ S ° Aa) ) 
\y/2S(mjJ ' 



W(5, S,5 0>c )d5 

(13) 



exp 



Then, 

dP(> to, a) 
dm 

' (14) 

and the Press-Schechter mass function can be found using 
equation (4), 



1 £q,c(q) dS 
~W 5(m) 3 / 2 dm 



S . c ( a ) 2 
' 2S*(to) 



dn 
dm 



(to, a) 



2 Pm,0 <5q,c(«) 



din y/S 



dim 



exp 



<^o,c(a) 2 
' 2S(m) 
(15) 



C. Derivation of the Double Distribution 

We now use the random walk formalism described in 
the previous section to derive the double distribution of 
the comoving number density of collapsed structures with 
respect to object mass to and local environment overden- 
sity Si, dn/(dmd8t). 

For the reasons described in the previous section, we 
replace the true overdensity field, S(x,a), with its linear 
extrapolation to the present time, 5(x,ao). Thus, we de- 
rive the double distribution of comoving n with respect to 
object mass to and extrapolated local environment over- 
density Si, dn/(dmdSi). We then use linear theory and 
the spherical evolution model to establish a conversion 
relation 5(5, a) and calculate dn/ '(dmdSi) as 



— — —(5i,m,a)dmd5i — — St(St,a),m,a dm—^dSt 

dmd5i dmdSi L - 1 odt 

(16) 

First of all, we need to define the local environment ex- 
trapolated overdensity Si in a precise way. We would like 
Si to be a measure of the density contrast of the medium 
in which a collapsed structure is embedded. Clearly, the 
value of 5i depends on how far from the structure itself 
its "environment" extends. We quantify this notion by 



introducing the clustering scale parameter, (3, which 
is defined in the following way: the "environment" of an 
object of mass to is a surrounding region in space which 
encompasses mass (3m (including the mass of the object). 
Hence, the local environment extrapolated overdensity 5e 
is the result of a filtering of 5(x, do) with a filter of scale 
(3m centered on the object. 

Formally, the above definition translates as follows. 
Consider the sharp in fc— space filtering function W m (k) 
discussed previously (eq. 6). The relation between the 
cutoff wavenumber and the filter mass, k c (m), is found 
by considering the form of the filter function in configu- 
ration space, 

_ sin[fc c (TO)r] - fc c (TO)rcos[fc c (TO)r] 

and multiplying by p m .o and integrating over all space, 
which yields [3] , 



fc c (m) 



67r 2 p m>0 



TO 



(18) 



For a collapsed structure at an epoch a which has mass 
to and is centered at a point x , we can write 

5(m,x ) = J W m (\x - x\) 5(x, a a )d 3 x = 5 , c (a) , (19) 

since the mean extrapolated overdensity of the collapsed 
structure itself is always the critical value for collapse, 
5 , c (a)- For that same object, the local environment ex- 
trapolated overdensity, 5(, is 

5 e (m,x ) = J W f3m (\x Q - x\)5(x,a a )d 3 x . (20) 

Equation (20) is then the definition of 5i for a given (3. 
In our double distribution, (3 is free parameter, which is 
however constrained to be between 1 and a few on phys- 
ical grounds. It cannot be < 1 since the mass of the 
object's environment always includes the mass of the ob- 
ject itself. In fact, as (3 approaches 1, the averaging which 
produces Si is taken only over the collapsed object itself, 
and inevitably returns the critical overdensity for col- 
lapse, 5o, c , for all objects. In the other extreme, /? > 1, 
the average 5 on a scale (3m is no longer a local quantity 
with respect to the central object. When (3 grows without 
bound, Si approaches for all collapsed structures, since 
averaging over the whole universe identically returns the 
background matter density, which corresponds to a van- 
ishing density contrast. In appendix A we show that 
our double distribution becomes proportional to a Dirac 
delta-function around St = in the limit (3 — > oo and 
proportional to a Dirac delta-function around St = <5o,c 
in the limit (3 — ► 1. 

We are now ready to use random walk theory results 
to calculate the double distribution. We first find the 
fraction of points in space which belong to structures of 
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mass between m and m+dm, which in turn are embedded 
in a medium of mean linearly extrapolated overdcnsity 
between Se and Se + dSe, f(m, Se, (3)d5e dm. The double 
distribution then is 



dn 



f(m,8e,f3)dmd5~e- 



dmdSe 

The quantity / can be written as 

fdmdSi = (f 1 dS e )(f 2 dm) 



(21) 



(22) 



where f\dSe is the fraction of points in space which have 
an average overdensity between 5g and Se + ddg on a 
smoothing scale (3m, and jidm is the fraction of points 
satisfying the previous condition which belong to col- 
lapsed structures of mass between m and m + dm. 

The first of the two factors above is the fraction of 
points still in the walk which are found between 5 and 
5 + dS at a "time" [3m. This is the solution of the ID 
random walk problem of 6g as a function of S, with an 
absorbing boundary at the critical collapse threshold <5o, c , 
as given by equation (12) but for a smoothing scale (3m, 



fidSe = 



exp 



' 2S(fim) 



— exp 



(Se-26 , c (a)y 



2S(/3m) 



^2irS{(3m) 



dS. (23) 



The second factor (^2) is the conditional probability 
that a point in space originating from (J3m, 5) in the mass 
- overdensity plane,will reach the "wall" for the first time 
for a smoothing scale between m and m + dm. This 
is then the probability that a particular point in space 
is absorbed by the "wall" So tC (a) at a particular "time" 
S(m), provided that the origin of the walk is transferred 
from (0,0) to (S(/3m),8e). This probablility can then 
be found if, in the expression for dP(> m,a)/dm (eq. 
14), we perform the substitutions So yC — > So yC — Se and 
S(m) — > S(m) — S(f3m). Similar conditional probabilities 
were originally calculated by [2] and [3] in the context of 
rates of mergers between halos. In our case, it is 



fidm 



5o,c( a ) 



exp 



(<5o,c(o)-<5«) 2 
'2[S(m)-S(/3m)] 



2tt [S{m) - S(f3m)} 



3/2 



dS 



dm 



dm . 
(24) 



Equations (21) and (22) then give 



dn 



dmdSg 



T (m,6e,(3,a) 



r 1 \ 1 ex P 
Pm,0 oo,c(a) - on 



\ $1 - 




2S(/3m) 


— exp 



(j g -2V c (a)) 2 
2S{fim) 



m 



2tt 



[S{f3m)} 1 / 2 [S(m) - S{[3m)f /2 



dS 



dm 



exp 



- \ 2 1 



' 2[S{m) - S(f3m)] 



(25) 



with 6 t < 6 ,c(a) and (3 > 1 so S(m) > S{f3m) [35]. 
Equation (25) is the double distribution we have sought 
and is the central result of this paper. Integrating 
dn/(dmd5e) over Se yields the Press-Schechter mass func- 
tion, as it should. The result is independent of the value 
of (3. We explicitly perform this integration in appendix 
A. 

Note that the functional form of our double distribu- 
tion is similar with that of the integrand used by [13] 
in their calculation of the cross-correlation between dark 
halos and mass using random walk theory, however the 
second variance of the field (corresponding to our S(f3m)) 
in their case refers to a fixed clustering radius and is in- 
dependent of object mass. 



D. Converting between Se and Se 

In appendices B and C we derive exact expressions for 
Si(St, a) in the case of the spherical evolution model, for 
an f2 m = 1 (appendix B) and an f2 m + Cl\ = 1 (ap- 
pendix C) universe (note however that all of the equa- 



tions we have presented up to this point are cosmology- 
independent, and can therefore be adapted for any cos- 
mological model). 

An excellent approximation to these conversion rela- 
tions can be derived from the expression 



5 a w 5 C 



l-a+^r 1 ^" . 



(26) 



Similar approximations were suggested by [23] and [15]. 
Equation (26) relates the linear overdensity at a time a to 
the true overdensity at the same time, and its accuracy 
is better than 2% throughout its domain for both £l m = 1 
and f2 m +f2A = 1 cosmologies. Its functional form is much 
simpler and more intuitive than the more accurate fit of 
[13]. The cosmological model enters only through S c . For 
the Einstein-deSitter universe, S c is given by equation 
(B17) for a co ii = 1, while for the f2 m + fl\ — 1 universe 
it is given by equation C29 for a = 1. Note that S c is 
related to the quantity So tC (a) (which appears explicitly 
in our double distribution expression) through 



(27) 
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where D(a) is the linear growth factor in the relevant 
cosmology. 

The limits of equation (26) are the same as the ones 
required for the exact conversion relation. When \8 a \ -C 
1, S a w 5 a . In addition, 5 a — > — oo as S a — » — 1, and 
<5 a — » oo as <5 — » Re- 
using equation (26), 



U(a) dc 



-l/«c 



(28) 



where a is the time at which we want to evaluate the dou- 
ble distribution. Note that close to virialization, equation 
(28) loses its applicability (as does the spherical collapse 
model), and has to be replaced by a recipe which does 
not diverge in 5. We present such recipes in appendices 
B and C, however the exact functional form of the con- 
version relation in this regime cannot affect any of the 
physically interesting results as the amplitude of the dou- 
ble distribution decreases rapidly enough with 5 that the 
contribution of the high-delta tail to the integrated mass 
function is negligible. We have verified this fact by com- 
paring the integral of our double distribution over 8 with 
the Press-Schechter mass function. When we extended 
the integration up to S c , the results agreed to the accu- 
racy of the numerical integration. When we extended our 
integration only up to S v (just below the application of 
our virialization recipe), the error relative to the Press- 
Schechter mass function was less than 0.02%. 



E. Clustering Scale Lengths and Correction for 
Central Object Contamination 

The definition of (3 and Se described above was suffi- 
cient for us to derive the double distribution from ran- 
dom walk theory. However, from a physical point of view, 
the presence of a collapsed structure at the center of the 
"environment sphere" contaminates the evaluation of the 
average "environmental" overdensity. If we want the dou- 
ble distribution to describe the properties of matter sur- 
rounding collapsed objects, we need to correct for the 
presence of the objects themselves. 

We will employ a simple, "top-hat" physical picture to 
calculate an appropriate correction (see also [15]). Note 
however that our correction is approximate, since the fil- 
ter wc used to smooth the overdensity field was k— sharp 
rather than top-hat in space. 

Let S c be the (true) overdensity of a collapsed object of 
mass m and radius i? v , and Se be the overdensity of the 
"environment sphere" of radius R c . The "environment 
sphere" encompasses a mass (3m, including the central 
collapsed object. We want to find the average overdensity 
S cxt of that part of the "environment sphere" which is 
external to the central object. For the collapsed object 
we can write 



to = -irRl{l + 6 c )p n 



(29) 



where p m is the mean matter density of the universe at 
the epoch of interest. For the environment sphere, in- 
cluding the central object, we can write 



(3m = -TrRl(l + 6e)p n 



(30) 



From equations (29) and (30) we get R% = R%(1 + 
5e)/(3(l + S c ). It thus follows that the length scale R e 
associated with the clustering parameter (3 is 



R Q 



P(l + S c ) \ 1/3 
1 + S e J 

3/?TO 

4tt(1 + 5 t )p 



1/3 



(31) 



We see that for a fixed clustering parameter (3, the 
length scale associated with an object of mass to is mass- 
dependent, scaling linearly with the virial radius but 
larger by a factor [(3(1 + 5 c )/(l + S e )] 1/3 > 1. Thus, wc 
can roughly think of the clustering scale parameter as a 
measure of how many virial radii we include as the local 
environment around each structure [36]. 

Having identified the environmental length scale, we 
can now isolate the environmental overdensity from that 
of the collapsed object. The volume of the environment 
sphere external to the central object contains a mass 



{(3 - l)m - -tt(R 3 c - Rl)(l + 8 cxt ) Pn 



(32) 



Using equation (31) to eliminate i? v , and dividing by 
equation (30) we obtain 



Q3-l)(l + fr)(l + g c ) 
<W ~ (3(l + 6 c )-(l + 6 e ) 



(33) 



which is the contamination-corrected overdensity for an 
environment sphere with uncorrected overdensity Se. 
Then, the contamination-corrected double distribution 
will be given by 



dn 



dmd5 C xt 
where 



((5 0Xt ,TO, a) 



dn dS e 
[o£(d ex t,a),m,a\ ■ 



dmdSg 



dd c 



t) = 



- 1 



and 



dSf 



(34) 
(35) 

(36) 



/?(l + gext)(l+<*c) 
(/?-l)(l + <y c ) + (l + *ext) 

= (3((3 ~ + S c ) 2 
d6 ext [(/3-l)(l + 5 c ) + (l + ,5 ext )] 2 

F. Derivative Quantities 



We now have enough tools to calculate derivative quan- 
tities of interest. The number density of collapsed objects 
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of mass greater than some minimum m m ; n [37] embed- 
ded in a medium of local overdensity between S ext and 
S cxt + d5 cxt is 



dn 



(> m min )d5 cxt = dS cxt 



dSe d5 c 



/>oo 
J m—r. 



dn 
a dmdSe 



-dm , 
(37) 

while the density of matter in collapsed objects of mass > 
m m i n embedded in a medium of local overdensity between 
S ex t an d ^oxt + dS cxt is 



dS, 



d P \jx jx dSe d6i 



L 



dn 

m —dm . 



dSe d5 cxt J m=mtain dmdS, 



(38) 

Of all the matter in the universe which belongs to col- 
lapsed objects of mass > m m i n , the fraction by mass 
which lives in underdense neighborhoods is 



rS 
J 6, 



-1 dS a: 



-(> m min )d5 cxt 



(39) 



Then, the mass fraction of the matter defined above 
which lives in overdensities will be /p, v = 1 — /p,un- 

Similarly of all the objects with mass m > m m ; n , a 
fraction by number which lives inside underdensities is 



/n, 



J<5cxt = — 1 



dn 
<i<5oxt 



(> m min )d6 e 



ft ext = -l di™ t ( > m min)cWext 



(40) 



The complementary number fraction of such structures 
living inside overdensities will be / niOV = 1 — /n.un- 

The number-density-weighted mean <5 0Xt for structures 
of mass > m m i n is 



^ = IsL=-l ^dS^O TO min)^ext 
Is^ t = - 1 dS^ ( > m min ) dS ext 

with a variance 



(41) 



2 



(42) 



Similarly the matter-density-weightcd mean S for 
structures of mass > m m ; n is 



(S) P = 



I 



-1 ^ oxt TO min)^e 



-1 dS a 



-(> m min )dS ext 



(43) 



with a variance 



dp 



(> m mi „)rf5ext 



III. RESULTS 

In this section we present plots of the double distribu- 
tion itself as well as of its various physically interesting 
derivative quantities. We compare results derived for a 
concordance, O m + Oa = 1 universe with WMAP pa- 
rameters (a s = 0.84, h = 0.71, n m = 0.27, il h = 0.04, 
[24]), and for an Einstein-deSitter (Q m = 1) universe with 
h = 0.71, n h = 0.04, but a 8 = 0.45. The different power- 
spectrum normalization in the Einstein-deSitter case was 
selected so that the Press-Schechter mass function in this 
case coincides with that of the concordance universe on a 
mass scale of 5.5 x 1O 14 M , which is between the values 
of m 8 (mass included in a sphere of comoving radius 8/i~ 4 
Mpc) for the two cosmologies (ms = 2x 10 14 M for the 
concordance universe while m 8 = 8 x 10 14 M Q for the 
Einstein-de Sitter universe). This value of as is also con- 
sistent with the fits of [9] given the WMAP result for the 
concordance universe. Finally, we use fitting formulae of 
[25] for the adiabatic cold dark matter transfer function 
to calculate the density field variance S(m). In this sec- 
tion, S always refers to <5 cxt , the true overdensity of that 
part of the "environment sphere" which is external to the 
central object. 

Figure 1 shows a 3-dimensional rendering of our double 
distribution as a function of mass and overdensity for 
fixed (3 = 2 and z = 0. The left panel corresponds to the 
concordance universe while the right panel corresponds 
to the Einstein-deSitter universe, and this arrangement 
is retained throughout this section. 

The features of the double distribution arc demon- 
strated in more quantitative detail in Figures 2-5. Fig- 
ures 2 and 3 show slices of the double distribution at 
constant values of mass. In Figure 2, different curves 
correspond to different values of the central object mass. 
In Fugure 3, all curves are for an object mass m = 
5.5 x 1O 15 M0, and different curves correspond to dif- 
ferent redshifts. Their most prominent feature is the 
pronounced peak at a relatively low value of \5\, indi- 
cating that for each given pair of z and to, there is a 
preferred, "most probable" value of the local environ- 
ment density contrast. As we can see in Figure 2, the 
location of this peak moves to higher values of the den- 
sity contrast as the mass of the object increases: small 
structures are preferentially located in relative isolation, 
while larger structures arc more likely to be found in clus- 
tered environments. This result fits well in the picture 
of hierarchical structure formation, as smaller structures 
tend to be merged into higher-mass objects as time pro- 
gresses. Lower-mass objects which are initially part of 
underdensities are less probable to undergo mergers, and 
hence are more likely to survive at late times than ob- 
jects which are initially part of overdensities. Conversely, 
higher-mass structures are more likely to be parts of over- 
densities where they can accumulate mass more easily 
through mergers with smaller structures. 

Note, however, that in the hierarchical structure for- 
mation picture, the mass scale where the exponential 
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FIG. 1: Surface plots of the double distribution for z — and {3 — 2 in fi m + £3a = 1 (left panel) and Einstein-deSitter 
(right panel) universes. The mass is measured in M . The vertical axis is linear, with the m — 5 axes level corresponding to 
dn/dmdSi = and the highest point corresponding to dn/dmdSe = 2.72 x 10 -2 (left panel) and dn/dmdSt — 1.89 x 10 _1 (right 
panel) objects per Mpc per 10 1j Mq. 



suppression of collapsed structures sets in increases with 
time. Thus, any given mass scale starts out as being 
a "high mass" at early times and eventually becomes a 
"lower mass" as it enters the power-law regime of the 
Press-Schechter mass function. Hence, according to the 
argument we used to explain Figure 2, the double distri- 
bution for any given mass scale should peak at increasing 
8 values with increasing redshift. This is because a partic- 
ular mass scale used to be closer to the high-mass end of 
the halo distribution in the past than it is today. Indeed, 
this is the trend seen in Figure 3. As we would expect, 
the peak of the distribution moves to higher 8 values with 
increasing redshift. The significantly more pronounced 
suppression of this mass scale in high redshifts in the 
Einstein-deSitter universe is due to the different power- 
spectrum normalization in the two cosmological models. 
Because of our choice in the power-spectrum normaliza- 
tion, the exponential suppression in the number density 
of structures sets in at lower masses in the Einstein-de 
Sitter case than in the concordance universe. Thus, there 
is a tendency to see more structures of higher mass in our 
concordance results than in the Einstein-de Sitter case, 
despite the intuitive expectation that a higher Q m uni- 
verse should have more massive structures at late times 
due to its ability to continue to form structures even at 
the present epoch. This would indeed have been the case 
if the power-spectrum had been normalized in the same 
way. 

That halos of a given mass are more strongly clustered 



with increasing redshift was also found by [26] , who used 
As (m) (the rms overdensity in the number of haloes more 
massive than some mass scale after smoothing with a 
spherical top- hat filter of comoving radius 8/i _1 Mpc) as 
a measure for halo clustering. A tendency of higher mass 
objects to be found in overdense regions was discussed by 
[13] and [12], who interpreted it by viewing halos today 
as progenitors of future larger-scale structures viewed at 
"high" or "low" redshift. 

In addition to the main peak at low \S\, an additional, 
much lower and sharper peak can be seen right before 
the critical overdensity cutoff. This peak is the result 
of the change of the functional form of the conversion 
relation between linearly extrapolated and true density 
contrast close to virialization, when application of the 
spherical collapse morel would lead 8 to diverge. The 
particular shape of the peak is an artifact of the recipe 
we adopted for dealing with the virialization regime, and 
carries no physical meaning (the shape of the peak is the 
shape of the high-<5 end of d5/dS). However, since the 
boundary conditions we use for 8(8) and its derivative 
are physical, we do expect to have some form of local 
maximum at the high-(5 end of the double distribution. 
Still, as discussed in the previous section, the effect of 
the details or even the existence of this local maximum 
on the physical quantities of interest is negligible. 

The high-<5 cutoff occurs at higher values of 8 in the 
concordance universe than in the Einstein-deSitter uni- 
verse. This is a result of the different density contrast 
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FIG. 3: Slices of the double distribution function at m = 5.5 x 1O 14 ,M and for various values of redshift z, for fi m + J1a = 1 
(left panel) and Einstein-deSitter (right panel) universes. The units of the double distribution are number of objects per Mpc 3 
per 10 15 M Q . 



achieved at virialization by structures in the two differ- 
ent cosmologies. In the Einstein-deSitter case this den- 
sity contrast is always 187T 2 , while in the concordance 
universe it is always higher and increases with time. At 
high redshifts, before the effect of A becomes significant, 
S c is very close to 187r 2 in the concordance universe as 
well, as can be seen in Figure 3. 

Figures 4 and 5 show slices of the double distribution 
at various fixed values of S, with z = and (3 — 2. Figure 
4 plots slices corresponding to relatively low values of |<5| 
(S = —0.5,0,0.5 and 3, close to the distribution peak in 



5). At the high-mass end of the distribution, the abun- 
dance of objects increases with increasing S, while in the 
low mass end of the distribution the trend is reversed, 
and the object abundance increases with decreasing S. 
This is in agreement with the behavior observed in the 
constant-m slices. 

Figure 5 plots slices corresponding to high values of 5 
(S = 10, 20 and 30), farther from the distribution peack. 
In this case, the curves do not cross, and an increase of 
5 simply results in an overall suppression of object abun- 
dance: structures of all masses are unlikely to be found 
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III 



FIG. 4: Slices of the double distribution function at constant values of 8 for z — 0, j3 — 2 and for fi m + Oa = 1 (left panel) and 
Einstein-deSitter (right panel) universes. Solid line: 8 = —0.5; dashed line: 5 = 0; dot-dashed line: 8 = 0.5; double-dot-dashed 
line: <5 = 3. The units of the double distribution are number of objects per Mpc 3 per 1O 15 M0. 




in 



FIG. 5: Slices of the double distribution function at constant values of 8 for z = 0, /3 = 2, and for O m + = 1 (left panel) 
and Einstein-deSitter (right panel) universes. Solid line: 8 = 10; dashed line: 8 = 20; dot-dashed line: 8 = 30. The units of the 
double distribution are number of objects per Mpc 3 per 1O 15 M0. 



overly clustered. This is because the final stages of col- 
lapse proceed rather quickly compared to the time spent 
around turnaround. The likelihood of a region observed 
in its late stages of collapse but before virialization is 
then low because the lifetime of this phase is small. 

In figure 6 we plot dn/dS(> 1O 12 M ) as a function of 
1 + S for different values of rcdshift. It is striking that at 
z = 0, the distribution peaks at negative 5 values (around 
5 = —0.6 in the concordance and —0.7 in the Einstein- 
deSitter universe), indicating that the most probable lo- 



cation for a collapsed object of mass > 10 12 M is an un- 
derdense environment. For the specific mass range, this 
trend is reversed by z = 1, when the preferred location of 
these objects is close to the universe mean (S = 0). This 
time-evolution pattern is independent of cosmology, as 
it is present both in the concordance and the Einstein- 
deSitter universes, and appears rather to be a charac- 
teristic of the hierarchical nature of structure formation. 
Parameters of this distribution can be calculated using 
equations (41) and (42) which, for the concordance cos- 
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FIG. 6: Distribution of structures with respect to local density contrast, dn/dS(> 10 1 M©), for /3 = 2 and for f2 m + Qa = 1 
(left panel) and Einstein-deSitter (right panel) universes. Solid line: z — 0; dashed line: z = 1; dot-dashed line: z — 2; 
double-dot-dashed line: 2 = 3. The units of dn/d5 are number of objects per Mpc 3 . 
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FIG. 7: Distribution of density of matter inside collapsed structures with respect to local density contrast, dp/dS(> 10 12 Mq), 
for f3 = 2 and for Q,m + Qa = 1 (left panel) and Einstein-deSitter (right panel) universes. Solid line: z = 0; dashed line: z — 1; 
dot-dashed line: z = 2; double-dot-dashed line: z = 3. The units of dp/dS are M per Mpc 3 . 



mology and z = give (<5) n = 0.43 and as, n = 4.36. 
The large value of the variance shows that the distribu- 
tion is significantly broad. However, the positive value of 
the mean is an artifact of the asymmetric boundaries of 
the distribution and its long high-<5 tail. This is demon- 
strated by the notably different locations of the mean 
and the median. The value of the latter is S = —0.22, 
therefore more structures in this range reside inside un- 
derdensities. 

Figure 7 is the matter-density counterpart of Figure 



6, as it plots dp/dS(> 10 12 M Q ) as a function of 1 + S 
for the same values of redshift as in Figure 6. Again, 
at the current cosmic epoch, the distribution peaks at 
negative values of S. Most of the virialized matter in the 
universe today appears to reside inside isolated objects 
rather than in clusters (note that decreasing the value of 
™ m in will only enhance this result since the trend towards 
isolation is more pronounced for the lower-mass objects). 
The trend of the peak with time (towards larger S for 
higher reshifts) is duplicated here as well. In particular, 
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FIG. 8: Fraction by number / n , U n (solid line) and my mass f P ,un (dashed line) of objects of mass > 10 12 Mq living in underdense 
regions, as a function of redshift, for O m + = 1 (left panel) and Einstein-deSitter (right panel) universes. 



note that at present, a significant fraction of the mass 
lies in moderately underdense regions. Equations (43) 
and (44) give for this distribution (in the concordance 
cosmology and for z = 0), (6) p = 1.20 and as lP = 6.23. 
The median of this distribution is at 5 = 0.20, a positive 
value. 

Finally, Figure 8 plots the evolution with redshift of 
the fractions by number and by mass, /n,un and /p ; uni of 
objects with m > 1O 12 M , living inside underdense re- 
gions. At high redshifts, when the mass of such objects is 
well above the exponential suppression cutoff, practically 
none of them is found inside undcrdcnsitics. This trend 
is reversed as the redshift decreases. In the ft + Oa = 1 
universe, an equal number of these structures are located 
inside underdensities by redshift 0.3 and by the current 
cosmic epoch, about 60% by number (but only 40% by 
mass) of these structures are located inside underdensi- 
ties. 

Figure 9 demonstrates the effect of changing the clus- 
tering scale parameter on the double distribution. Slices 
of the double distribution along m=5.5x 1O 14 M are 
plotted (in linear axes) as a function of S, and for (3 = 1.5 
(solid line), 2 (dashed line), 3 (dot-dashed line) and 10 
(double-dot-dashed line). The location of the peak ap- 
pears to be extremely insensitive to the value of (3 for 
moderately low values. It very slowly moves towards 
S = with increasing /?, as it should (increasing [3 re- 
sults in averaging the overdensity over increasingly large 
volumes). Note that even as (3 approaches 1, the peak 
will not move towards S c , as a result of our correction for 
the central-object contamination. This makes our for- 
malism particularly suitable to study the properties of 
matter very close but outside a virialized structure (e.g. 
the local density of accreted gas). 

The effect of j3 on an integral quantity is shown in 



Figure 10, which plots dn/dS{> 10 12 M Q ) for (3 = 1.5 
(solid line), f3 = 2 (dashed line) and (3 = 10 (dot-dashed 
line). Again, the results are extremely insensitive to the 
value of (3, which gives us confidence about the robustness 
of the location of the peak of our distributions. 



IV. DISCUSSION 

We have presented an extension of the Press-Schechter 
mass function, in the form of double distribution of struc- 
tures with respect to mass and local overdensity. We 
have done so by introducing a clustering scale param- 
eter (3 > 1, which we use to associate with each col- 
lapsed object of mass m a larger environment of mass 
(3m. The scale parameter (3 can be expressed as a fuc- 
ntion of the number of virial radii included in the local 
environment of each structure. We found that for rea- 
sonable values (3 ~ 2, the shape of the distribution does 
not depend sensitively on this parameter. Integration 
over linearly extrapolated overdensity returns the origi- 
nal Press-Schechter mass function, independently of the 
value of (3. 

We present the double distribution in terms of the true, 
physical, nonlinear density contrast 5. However, in cal- 
culating the distribution it is useful to identify regions 
using instead the overdensity obtained via linear analy- 
sis, 5, extrapolated to the present epoch. A useful fitting 
function is given for the 6 — S conversion. 

The double distribution is useful because it allows us to 
have an explicit analytical if approximate description of 
the environment in which collapsed objects of all masses 
reside. Using the tools we have developed, it can be 
readily calculated for any flat cosmology, and evaluated 
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FIG. 9: Slices of the double distribution function at m = 5.5 x 10 14 , Mq and for different values of the clustering scale parameter 
j3, for 57 m + Qa = 1 (left panel) and Einstein-deSitter (right panel) universes, plotted in linear scale. Solid line: (3 = 1.5; dashed 
line: (3 — 2; dot-dashed line: (3 — 3; double-dot-dashed line: (3 = 10. The units of the double distribution are number of objects 
per Mpc 3 per 1O 15 M . 
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FIG. 10: Distribution of structures of mass larger than 10 12 Mq with respect to local density contrast, dn/dS, for fi m + Q,a = 1 
(left panel) and Einstein-deSitter (right panel) universes, at z = 0, and for f3 = 1.5(solid line), (3 = 2 (dashed line) and 
(3 = 10(dot-dashed line). The units of dn/dS are number of objects per Mpc 3 . 



at any epoch. Consequently, it offers new insight into the 
growth of structure as well as the present distribution of 
collapsed objects. 

We have evaluated the double distribution and some 
of its integral moments for both a concordance cosmol- 
ogy and an Einstein-de Sitter universe. Some key results 
are that at any redshift, the double distribution is dom- 
inated by a peak which shifts in mass but is always at a 
relatively low value of \S\. For each mass, there is a most 



probable 6, which increases with structure mass. More- 
over, at the present epoch in the concordance universe, 
the most probable environment is a modest wnderdensity, 
for all objects below about 1O 14 M ; thus, underden- 
sities are preferentially populated by low-mass objects. 
Finally, the fraction of mass in undcrdensities increases 
with time, and in the concordance cosmology the present 
underdense mass fraction in objects of M > 10 12 M Q is 
about 40%. These trends can be understood in terms of 
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hierarchical clustering in which overdense regions are the 
site of vigorous merging that clears out low-mass objects, 
which then find their last refuge in voids. 

These results are consistent with other analyses in 
the literature which use Press-Schechter-like formalism 
to probe the correlations between structures. For exam- 
ple, [13] extend the Press-Schechter random- walk picture 
in a way very similar to ours, but identify a local envi- 
ronment in terms of a fixed, mass-independent radius. 
Our results are thus complementary to these, because we 
adopt a mass-dependent environment based on the virial 
radius. 

Finally, it is well-known that the Press-Schechter mass 
function provides an excellent characterization of the re- 
sults of numerical TV-body simulations (e.g., [8]). It 



would be of great interest to compare the analytic dou- 
ble distribution we have presented with numerical results. 
We plan to address this issue in future work. 
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lower bound, not only due to the finite mass of the dark 
matter particle, but also due to the existence of a dark 
matter Jeans mass, however small this may be. Practi- 
cally, the Press-Schechter dn/dm diverges asm^O and 
setting a minimum mass is required to extract interesting 
information. 
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APPENDIX A: ANALYTICAL PROPERTIES OF 
THE DOUBLE DISTRIBUTION 

1. Derivation of the Press-Schechter Mass Function 
From the Double Distribution 

Using Si to denote S(m) and 5*2 for S(fim) we have: 
I 



5o - - dn 

do 

oo dmdS 



Pm 

m 



dSi 



dm 



1 



p ' 




(<5o,c 
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2S 2 
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-S 2 ) 
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m 2ttSI /2 (Si- S 2 ) 3 / 2 
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7T TO 2 y/Sl 
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S dS exp 



s± — s 2 
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2{Si - 5 2 ) 



25 2 



exp 



S'2 



2(5! - 5 2 ) 



Si 



3/2 



y r 2^S 2 exp 



25i 



(Al) 



where we performed a change of variables 8' = 8o, c — 8, 
and we set 8' — > —5' in the second integral. The final 
result is the Press-Schechter mass function formula. Note 
that this result is independent of the value of (3. 



Behavior of the Double Distribution in the limit 

(3 -f oo 



Then, using the notation of the previous section, 



lim dn r (m,S e ,(3,a) = lim dn (Si, S 2 , 8 e , a, m) . 
P^oc dmdde i> 2 ^o drnddg 

(A3) 

The limit of a unit-area Gaussian when its width vanishes 
is the Dirac delta- finction Sr>, 



In order to find the behavior the double distribution as 
P — > oo, we recall that, because 5(m) decreases mono- 
tonically with m, its limit in the infinite f3 regime will 
be 



lim 



cxp 



2A 2 



= 6u(x-x ). (A4) 



lim S((3m) = . 



(A2) Using this result, we get 



lim 



dn 



S2->o dmdSi 



Pm 

m 

Pm 

m 

Pm 

m 



dSi 
dm 

dSi 
dm 

dSi 
dm 



8o,, 



2ir 

kc - Se exp 



lim ^ exp 

S2 — ^0 



(*o,c - kf 



2(5i - 5 2 ) 



exp 







2S 2 _ 


— cxp 



($e — 2<5 0iC ) 2 



2S 2 



2Si 



'2tt 5 x ° 
*o.c - ^ ex P f 



3/2 



lim 





*< 


exp 


2S 2 



5 2 1/2 (5i -5 2 ) 3 / 2 
exp 



— lim 

s 2 ^o 



(^-2^ . c ) 2 
2S 2 



2Si 



2tt 



5r 



3/2 



s 2 ^o V2?r5 2 
£d(<^) - 5d(^ - 25 ,c) 



V2T5^ 



(A5) 



16 



However, the 5p— domain of the double distribution is 
between — oo and <5o, c , and therefore the value 5p = 2<5o, c 
is outside its domain. Hence the second Dirac delta- 
function is always zero, and 



dn 

hm — 

f>^°° dmdSi 



Pm_ 

m 



dm 



£>o,c 



cxp 



(S ,c-S e f 
2Si 



2tt 



s 3/2 



Sd(Si), 
(A6) 



proportional, as expected, to a Dirac delta-function cen- 
tered at Se = 0. 



3. Behavior of the Double Distribution in the limit 

15-1 

Denoting S((3m) by S2 and letting <fi = S(m)/S([3m) 1 
we seek the behavior of the double distribution in the 
limit (3 —> 1 or <\> — > 1. Dchning 



(the linear-theory density contrast at epoch a) should 
agree at early times. 



Spherical Evolution Model in an fi n 
Universe 



The evolution of a spherically symmetric, overdense 
perturbation in an otherwise homogeneous Q m = 1 uni- 
verse is described by the parametric equations 



Pm So, c — 
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we can write 
dn 



lim 



*! dmddi 
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(1 — cos 9) , and a = a co n 



— sm( 
2tt 



2/3 



(Bl) 

where a co ii is the scale factor of the universe when the 
perturbation formally collapses to a point, a p is the scale 
factor of the perturbation, and 9 is the development an- 
gle. Note that the perturbation will turn around (reach 
its maximum size, a Pimax = 4a co n(127r)~ 2 / 3 ) when 9 = n, 
at a time a = a co n/2 2 / 3 . 

The normalization of equation Bl is such that the den- 
sity contrast 5 can be expressed as 



5=1- 



1 . 



(B2) 



Hence, for any density contrast 5, equations Bl and B2 
can be combined to give a unique development angle 6(6) 
which is the solution to the transcendental equation 



6 2 / 3 (fl-sinfl) 2 / 3 
2(1 - cos 6>) 



+ 



1/3 



0. 



(B3) 



(6o,c - 5 e y 



proportional, as expected, to a Dirac delta-function 
around 6q c . 



Similarly, the spherical evolution solution for an un- 
dcrdcnsity is given by the parametric equations 

g2/3 

a p = A p (cosh 77 — 1) , and a = A p — — (sinhry — rj) 2 ^ 3 . 

(A8) " (B4) 

where r\ is the development angle in this case. Equation 
B4 together with equation B2 can be combined as before 
to give r\(6) as the solution to the transcendental equation 



APPENDIX B: RELATING LINEARLY 
EXTRAPOLATED AND TRUE OVERDENSITIES 
IN AN EINSTEIN-DE SITTER (fi M = 1) 
COSMOLOGY 

In this appendix we derive a conversion relation 6o(a, 5) 
for an f2 m = 1 cosmology (here, 6 is the density contrast 
predicted for a density perturbation at cosmic epoch a by 
the spherical evolution model and 6q is the extrapolation 
of the density contrast to the present cosmic epoch using 
linear theory). In order to do so, we hrst calculate 5(a) 
from the spherical evolution solution, then calculate So 
using linear theory, and finally require that 5(a) and 5 a 



2. S (a, 8) according to the spherical evolution 
model 

The behavior of 5 in the linear regime in this cosmology 



is 



5 = 5ga . 



(B6) 



This result should coincide with the linear expansion of 
the spherical evolution result at early times. Expanding 
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the parametric solution to second nonvanishing order in 
and eliminate 0, we obtain 



a p (a) = a 



(12tt) 2 / 3 a 



20 



dcoll 



(B7) 



We then substitute equation B7 in the definition of 6 (eq. 
B2) to get 



j 3(12tt) 2 / 3 

o = — — a 

20a co ii 

which, by comparison to equation B6 gives 

5 3(12tt) 2 / 3 
°° = — ™ • 

Then, the conversion relation we seek is 

6 2 / 3 3 



6 (a,S) 



20a 



[9(5) - sin 0(6)] 



2/3 



(B8) 



(B9) 



(BIO) 



where 0(6) is given by equation B3. 

Equation BIO has the undesirable property that it di- 
verges as approaches 2ir. This is of course a conse- 
quence of the perturbation formally collapsing to a singu- 
larity in the spherical evolution model instead of reaching 
virial equilibrium. If wc make the usual assumption that 
at virialization the radius of the perturbation is a max /2 
and we additionally require that 

• S (a,S) is continuous and smooth at = 3n/2 

• a p — a p , max for all a > a coU 

then for > 3ir/2 (which corresponds to 6 > 9(3ir + 
2) 2 /8) we can replace equation BIO with 

6 (a,6) = 6 , v + 6' v (6 - 6 V ) 

3(5o, c -5 , v )-(5 c -6 v )(26i v + 6i c ) 



+ 



(S C - Syf 



(6 - 5 V ) 2 



(K, c + fyv)(<?c - S v ) - 2gp^ - SqJ 



(Sc - Sy) 3 



(Bll) 



(see appendix C for a discussion of the reasons for em- 
ploying this particular functional form, and §11 and III 
for a discussion on why the effect of such a choice on the 
double distribution is negligible). In equation Bll, 
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10o(l + <5 C ) 2 / 3 



(B12) 



The last equality coming from the fact that after a co ii 
the radius of a perturbation remains constant and equal 
to a Pimax /2, while its density contrast 6 changes only 
due to the expansion of the background universe, S — 
(2a/a Pimax ) 3 - 1 or S = (10a<5 /3) 3 - 1. Finally, 6 , c is 
given by equation B17 while S' v is given by equation B18 
for 6 = 5 V and — '3tt/2. 

To get the linear behavior of 5 for an undcrdensity we 
expand the parametric solution B4 to second nonvanish- 
ing order in r] and we eliminate r] to get 



6 2/3 



1 + 



1 a 
10 X 



(B13) 



Substituting equation B13 in the definition of S (eq. B2), 
we get for the time dependence of 6 at early times, 



S = 



10A C 



(B14) 



from which, by comparison to equation B6, we get 

s 3 



'Jo 



10A n 



Then, So(a,5) will be 



5 (a,S) 



6 2 / 3 3 
20a 



[sinhry((5) — r](S)] 



2/3 



(B15) 



(B16) 



where rj(5) is given by equation B5. 

Equation B16 is valid for all rj and its limit as 5 — > — oo 
is 6(6) — > — 1. Thus, although the linearly extrapolated 
field can become < — 1, the corresponding value of the 
actual 6 is always > — 1, as the physical requirement 
P P > demands. 



3. Critical extrapolated overdensity for collapse, 

So, c (a) 

The critical extrapolated overdensity for collapse can 
be found from equation B9 



<5o,c(ocoii) — 



3(12tt) 
20" 



2/3 



-a col] rj 1.69a coll . 



(B17) 



Note that the above equation has the functional form 
<Vc(acoii) l/-C( a coii), where D(a) is the linear growth 
factor for this cosmology. This is also true in the f2 m + 
Cl\ = 1 case. 



4. ddo/dS\ a 

In addition to the relation between 6 and So, we will 
also need the derivative d5o/dS\ a in order to convert be- 
tween true and extrapolated overdensity differentials in 
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equation (16). In the case of an overdense structure, 
S > 0, equation BIO gives 



dSo 
05 



6 2 / 3 l-cos6»((5) dO 



(B18) 



10a [6(d) -sin 0(5)] 1/3 dS ' 

To evaluate dO/dS we define the auxiliary function 

F a (6,5) = 6 2/3 (0 - sintf) 2 / 3 - 2(1 - cos0)(l + <5) 1/3 . 

(B19) 

From equation B3 we get immediately F a (9, 5) = 0, and 
differentiating we get dF a = = ^§fd0 



^ffdS. Hence, 
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Equation B18 is valid only for < 5 < 5 V . For 5 > S v 
equation Bll gives 
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In the case of an underdense structure, S < 0, equation 
B16 gives 



(Mo 
d5 



6 2 / 3 cosh7?(£)-l dr] 
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As before, in order to evaluate dr]/d5 we define the aux- 
iliary function 

G a (?7, 5) = 6 2/3 (sinh77 - t]) 2/3 - 2(coshr/ - 1)(1 + S) 1/3 . 

(B25) 

Equation B5 implies G a (r],5) = so 



where 



and 
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APPENDIX C: RELATING LINEARLY 
EXTRAPOLATED AND TRUE OVERDENSITIES 
IN AN SIm. + fi A = 1 COSMOLOGY 

In this appendix we will derive a conversion between 
true and extrapolated overdensity, 5 (So, a) for an £! m + 
Cl\ = 1 cosmological model. We will do so by first cal- 
culating the true density contrast 5(a) of a density per- 
turbation at cosmic epoch a as predicted by the spher- 
ical evolution model, then calculating <5 , which is the 
overdensity of the same spherical perturbation if extrap- 
olated according to the linear theory until the present 
cosmic epoch, and finally requiring that at early times 
linear theory and the linear expansion of the spherical 
evolution model should give the same result. 



1. Spherical Evolution Model in an fi m + Qa = 1 
Cosmology 

a. The Evolution Equation 

In the spherical evolution model, the spherical density 
perturbation under consideration behaves as an indepen- 
dent non-flat sub-universe. Its evolution is dictated by a 
Friedmann equation, 



(^) 2 = ^o 2 ^a 2 (a- ; 



(CI) 



where a p is the radius of such a spherical density per- 
turbation in an otherwise homogeneous universe, u> = 
f^A/^m = — 1 (where f2 m and fl\ are the matter 
and vacuum density parameters of the background uni- 
verse) and k is a constant characteristic of the amplitude 
and sign of the perturbation: the larger the the larger 
the deviation from homogeneity at a given time, while a 
positive k corresponds to an overdensity and a negative k 
to an undcrdensity. Clearly then in equation CI, the first 
term in parentheses on the RHS is the matter term, the 
second is the vacuum term and the third is the curvature 
term, which can have a positive or negative sign depend- 
ing on whether we are studying an "open" (undcrdensity) 
or "closed" (overdensity) perturbation. The normaliza- 
tion of a p is such that, had the specific spherical region 
begun its evolution with no curvature (k = 0), a p at the 
present cosmic epoch would have been a p i K —, 



p(k=o),o 



1. For 



this reason, the density contrast 5 of the perturbation at 
epoch a is given by equation B2 

The behavior of the perturbation radius func- 
tion of the universe scale factor a can be found by taking 
the ratio of the Friedmann equations of the perturbation 
and the background universe, thus obtaining [27] 



da p 
da 



a Lua* 



+ Loa 2 



+ (C2) 



Lua 3 + 1 



Equation C2 implies that the smallest positive perturba- 
tion which will turn around and collapse corresponds to 
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the smallest positive k for which the equation 



uoa" 



Kfflp + 1 = 



has a real positive solution [9]. This gives 

Kmin.coll = 3cj 1 / 3 /2 2/3 . 

Equation C2 can then be re-written as 

1/2 



(C3) 



(C4) 



da p 
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/ -1 i 2 \ J-/- 2 



K < K m in,coll O r 

K _ K min, colli a < a ta 



1/2 



K _ K min. colli a > Ota 

where ata is the scale factor of the universe when the per- 
turbation reaches its maximum (or turnaround) radius. 
The turnaround radius is the smallest of the two positive 
solutions of equation C3, 





-3 
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\0J3 J 







(C6) 

Equation C6 has an asymptotic behavior a Pita ~ 1/k 
when k/cj 1 / 3 » 1, as expected from equation C3. 
The maximum possible turnaround radius, o P! ta,max is 



achieved for n 



^min.coll 



and is a 



p,ta,max 



(2a;)- 1 / 3 . All 



other collapsing overdensities will have a P:ta < a p ,ta,i 



b. Qualitative Description of the Evolution of Structures 

The introduction of the additional vacuum term in the 
Fricdmann equation considerably complicates the simple 
classification of density perturbations to overdensities (all 
of which turn around and collapse in an f2 m = 1 cosmol- 
ogy) and underdensities (all of which expand forever). 
In the fi m + J1a = 1 universe there exist overdensities 
which will continue to expand forever. The behavior of a 
perturbation in such a cosmology is parametrized by the 
quantity k/cj 1 / 3 , and we can identify the following cases. 

Case I, k/u} 1 / 3 < —1: large underdensities, ex- 
panding forever. The table below shows the relative 
magnitude of the three terms in the Friedmann equation 
(matter, curvature and vacuum term) for different val- 
ues of the scale factor of the perturbation. The first line 
in the table indicates the hierarchy of the three terms, 
from largest to smaller, for each range of the scale factor. 
The second line indicates the dominant term in each scale 
factor range. The third line shows the approximate de- 
pendence of the radius of the perturbation, a p , on time, 
assuming that only the dominant term contributes to the 
Friedmann equation in each range. 



°P < R 


R < a P < ife 


-J— < a < -i/i^i 


"P ■> V w 


MCV 


CMV 


CVM 


VCM 


matter 


curvature 


curvature 


vacuum 


a p ~ ei* 


a p ~ t 


a p ~ t 


a p ~ e* 



Case II, —1 < k/uj 1 / 3 < 1: small perturbations, 
expanding forever. These can be either underdensities 
(k < 0) or overdensities (k > 0). In both cases the 
curvature term never becomes dominant. The following 
table shows their different evolutionary stages (as in Case 
I)- 





\/^<«P<lfe 


W < ° p < l«l 


«p > R 


MCV 


MVC 


VMC 


VCM 


matter 


matter 


vacuum 


vacuum 


a P ~ t 2/3 


a p ~ i 2 / 3 


a p ~ e* 


a p — e* 



Case III, 1 < k/uj 1 / 3 < 3/2 2 / 3 : "coasting" 
overdensities, expanding forever. These overdensi- 
ties continue to expand forever despite the fact that they 
go through a phase in their evolution when the curva- 
ture term becomes dominant and their expansion slows 
down. During this phase, the contributions of the matter 
and vacuum terms, which are the ones driving the expan- 
sion, add up to a value always higher than the curvature 
term, although the curvature term is larger than each one 
of them. When the perturbation enters the curvature- 
dominated phase, the expansion rate decreases and the 
perturbation grows much more mildly than t 2 / 3 . The ex- 
pansion rate reaches a minimum at a p = (2w) -1 / 3 , after 
which it increases again as the perturbation approaches 
the phase of exponential expansion. This phase between 
the matter-like expansion and the exponential expansion 
is denoted by (*) in the table below. 



« P < i 




is < «p < Vz; 


«p> 


MCV 


CMV 


CVM 


VCM 


matter 


curvature 


curvature 


vacuum 


a p ~ 


(*) 


(*) 


a p ~ e* 



Cases I-III are all sub-cases of the Lemaitre model 
([28], [29]), which features an inflection point at a p . c = 
(2cj) -1 / 3 where d p = while d p > 0. The rate of expan- 
sion initially decreases to achieve its minimum (positive) 
value when a p = a p , c , after which point the expansion 
accelerates again. 

Special Case, k/uj 1 / 3 = 3/2 2 / 3 : Eddington 
Overdensity. This overdensity is the lowest k over- 
density which does not expand to an infinite radius. 
However, it does not turn around and collapse, but 
it approaches its (finite) turnaround radius, a p , max = 
(2U))- 1 / 3 (from eq. C6) as t -> oo. As seen by an observer 
inside this overdensity, as t — > oo the part of the universe 
outside a Pimax will accelerate away and eventually exit 
the horizon, and the observable universe ( "local Edding- 
ton bubble") will asymptotically approach the Einstein 
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TABLE I: Characteristic times of the spherical evolution 
model in an r^ m -t- Oy\ — 1 universe. 



event 

turnaround 



time 

tt 



matter-vacuum 
equality 



tMV 



fo 13 d % 



3^ — Kflp + l 



matter-curvature 
equality 



curvature- vacuum 
equality 



static universe (as in the Eddington model with a cosmo- 
logical constant). 

Case IV, k/u> 1/3 > 3/2 2/3 : large overdensities, 
eventually collapsing. When such a structure enters 
the dominant-curvature-term phase, its expansion rate 
starts to decrease (a p ~ t e with e = e(t) monotonically 
decreasing from 2/3 to 0), until the expansion halts, at 
a p = a Pjta which occurs at a time ita, given in table I. 
After t ta the perturbation turns around and contracts, 
its evolution being symmetrical in time about t ta , i.e. 
a p (t) = a p (2tt a ,—t) for t > tt& (this is a consequence of eq. 
CI and holds for any cosmological model as long as the 
RHS of the Fricdmann equation involves no explicit time- 
dependence). Eventually, the perturbation will formally 
collapse to a singularity at time t co u — 2t ta . 
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CMV 


CMV 


MCV 


matter 


curvature 


curvature 


matter 




ftp ~ t e 


a P ~ (2tta- t) £ 


a p ~(2t ta - t]^ 


expansion 


expansion 


contraction 


contraction 



In all of the cases discussed above, the transitions be- 
tween different phases of their evolution occur at char- 
acteristic times, those of matter- vacuum equality tMV, 
matter-curvature equality £mc an d curvature-vacuum 
equality t C y. At these times (shown in table I), the 
corresponding terms in the Friedmann equation become 
equally important. Note that in the case of the Edding- 
ton overdensity and of case IV collapsing overdensities, 
matter-vacuum equality and curvature-vacuum equality 
are never reached, and the vacuum term never dominates 
over any of the other terms. 

In the next section we derive exact solutions for the 
time-evolution of a p for perturbations of different curva- 
ture. However, surprisingly accurate approximate solu- 
tions can be derived using only linear theory and equation 
(26). Solving for the spherical collapse density contrast 
we get 
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Since o p = a(l + 5 a ) 1 / 3 , we can write for collapsing 
overdensities 
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where the initial conditions (curvature) of the perturba- 
tion are parametrized by its collapse epoch, o co ii, while 
the cosmology enters through the functional form of the 
linear growth factor and the linear collapse overdensity, 
S c . Similarly for perturbations which expand for ever we 
can write 



S D(a)/D(a ) 



-i <5 c /3 



S D(a) 
5 C D(a ) 



S c /3 



(C9) 



where the curvature of the perturbation is parametrized 
by its extrapolated linear density contrast at the present 
epoch, So- Note that for overdensities which expand for- 
ever, So > and a p < a, while for underdensities Sq < 
and a p > a. Also, because D(a) asymptotes to a constant 
value for a — > oo (as we will see in the next sections) , 
a p grows proportionally to a at late times. This is the 
exponential expansion phase, described in our analysis 
above. 



c. Solutions of the Evolution Equation 

For eventually collapsing structures (k > K m in,coii)> 
separation of variables in equation C5 and integration 
yields, 
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(CIO) 

where at a is the cosmic epoch when a p = a p . ta - Now the 
integral on the LHS of equation CIO can be calculated 
using [9] 
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' 1/2 sinh" 



vW + 1 3 

The integral of the RHS can be re-written as 

2 



(Cll) 



/ 

Jo 



^fxdx 



o V ujx 3 — kx + 1 3 



^Vifcj*) (C12) 



where Vi is the incomplete vacuum integral of the first 
kind, defined in appendix D, and 



Op/ <^p,ta > 



M = (w4,ta) 1 • 



(C13) 
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Note that for this range of curvature values, k/uj 1 / 3 
(which is the quantity which parametrizes the behavior 
of the perturbation with time) is a function of n alone, 
with 



k/u 1 ' 3 = {l + l l)/ l l 2 l 3 . 



(C14) 



Using equations Cll and C12, equation CIO can be 
rewritten as, 



a < a ta 



a = 



w- 1 /3{ s i I1 h[2Vi(l,A*) - Vi(r, M )]} 2/3 , a>a ta 

(C15) 

Equation C15 is the spherical evolution solution for a 
collapsing perturbation with amplitude k (parametrized 
above by /i) in an fl\ + Q m = 1, fl\/fl m = u> universe, 
and gives r (and hence a p ) as a function of a for this 
model. Note that Vi(r,n) is the development angle for 
this cosmology. 

The scale factor of the universe at turnaround for a 
given collapsing overdensity can be found immediately 
from equation C15, 



at^tu- 1 / 3 [sinhVi(l, M )] 2/3 • 



(C16) 



The scale factor of the universe at collapse, a co ii (when 
the scale factor of the perturbation becomes formally 
zero, a p c = r c = 0) is, from equation C15 and since 
Vi(0, M ) = 0, 



a coll 



w- 1 / 3 [sinh2V 1 (l,M)] 2/3 



(C17) 



Equation C15 should not be applied "literally" until 
the final collapse of the perturbation to a singularity, 
since the physical picture for the late stages of the evo- 
lution of a perturbation involves virialization at a finite 
radius. It has been shown by [30] that the analogous 
arguments which give a PtV — a p , ta j2 for the f2 m = 1 cos- 
mology give, for an Q m + Cl\ = 1 universe, 



Auoa 



2 + 2uja 3 n 



'p.ta 



a p ,v + 1 = 0. 



(C18) 



The physically meaningful solution of C18 which gives 
the correct behavior for uo — > is (using eq. C13) 



a p,v — a p,ta 



2)1 



2 1 
- cos - 
3 3 



cos 



27^ 2 



(2/x + 2)3 



(C19) 



The scale factor a v of the universe when the scale factor 
of the perturbation past its turnaround becomes equal to 
will be given by the second branch of equation C15, 



*p,v; 



for a p = a Pl v Then, the validity range for the second 
branch of equation C15 is a ta < a < a v . 

For a > a v , we can no longer use the spherical evolu- 
tion solution to describe the physical picture of interest 



(virialization). In the next section we will present a sim- 
ple recipe we will use to follow the late stages of evolution 
of the perturbation which satisfies the desired boundary 
conditions (a p = a P;V at a co n and constant thereafter). 

For perpetually expanding structures (n < K m in,coii), 
equation C5 gives 



Jo J ujv 3 + 1 Jo 



y/xdx 



(C20) 



y/ Ujy 3 + 1 JO V W 3 - KX + 1 

In this case, the integral on the RHS can be rewritten as 
#= = f 1/2 «i(^) ( C21 ) 

° — KX + 1 O 



UJX' 



where Hi is the hyperbolic vacuum integral of the first 
d, defined in appendix D, and 



a P /|a p ,R.| , 



w = Ha p , R | 3 ) 1 



(C22) 



where a Pi R is the only real (and always negative) root of 
equation C3 when k < ft m in.coii, 



Op.R 



-U)3 



(C23) 

As in the case of collapsing perturbations, k/uj 1 / 3, is a 
function of vj alone, with 



k/lo 



1/3 



(1 — w)/w 



2/3 



(C24) 



Then, w = 1 is a flat subuniverse (not perturbed with 
respect to the background), and perturbations with vj > 
1 are underdensities while 1/4 < vo < 1 correspond to 
non-collapsing overdensities. 
Then, equation C20 becomes 



-1/3 



{sinh [Hi{r,zv)]} 



2/3 



(C25) 



which is the spherical evolution solution for a non- 
collapsing perturbation with amplitude k (parametrized 
above by w) in an Q\ + Q m = 1, Q\/fl m ~ u> universe. 
As for collapsing perturbations, Hi(r,zu) is the develop- 
ment angle. 



2. 8 (S) according to the spherical evolution model 

The linear theory result for a growing-mode perturba- 
tion in an S! m + Oa = 1 cosmology is [31] 



s=~s D{a) 



'D(a ) 



(C26) 



where D, the linear growth factor, is given by D(a) 
A[(2w) 1 / 3 a] with 



A{x) 



(x 3 + 2)V2 r x 



x 3 / 2 J \u 3 + 2 



3/2 



du . (C27) 



22 



To find the relation between 5o and k, we expand the ex- 
act (6 = a 3 /ftp — 1) and linear relations for the overden- 
sity to first order in a and demand that the coefficients 
be equal. Thus we get [9] 



(2w) 



1/3 



8 



(2w) 



1/3 



3A [(2w)V3ao] 3A [(2w)V3 



-<5 • 



(C28) 



since ao = 1 . This is the value of the constant k for 
a perturbation which at a cosmic epoch ao — 1 has a 
linearly extrapolated overdensity Sq . Note that since the 
linear theory result is the same for both underdensitics 
and overdensities, equation C28 holds for both cases. For 
underdensitics, both k and S will be negative, while for 
overdensities both will be positive. 

At any given epoch a, there is a unique perturbation 
(parametrized by k or, cquivalently, by /i or w) which 
will have achieved a true density contrast 5 at that time. 
Therefore, to calculate the desired conversion relation 
5o(a, S) we first calculate k (or, equivalently, fi or w) 
from the given a and 6 and the appropriate solution of 
the evolution equation (CI 5 or C25). Then, we use equa- 
tion C28 to evaluate 8q. 

To determine which is the appropriate solution of the 
evolution equation we need to use for each 5, we calculate 
the limits of applicability of each equation in terms of 6. 

• Equation C25 is applicable for all forever expand- 
ing perturbations. For any given epoch a, the 
maximum density contrast of such perturbations 
is achieved by the Eddington perturbation and is 
equal to fed (a), given in line 1 of Table II. Then, 
the applicability domain of equation C25 is — 1 < 
8 < fed (a) , and the conversion relation in this case 
takes the form shown in the 1st line of Table III. 

• The first branch of equation C15 is applicable 
for eventually collapsing perturbations which, how- 
ever, have not reached their turnaround radius yet. 
The maximum 5 of all such perturbations at a given 
a is achieved by the perturbation which is turning 
around at a, and is equal to <5ta(a), given in line 2 
of Table II. The applicability domain of the first 
branch of equation C15 is then fed (a) < 5 < <5 ta (a) 
and the conversion relation in this case is shown in 
the 2nd line of table III. 

• The second branch of equation C15 is applicable for 
eventually collapsing perturbations which are past 
their turnaround but which have not yet reached 
their virial radius. The maximum S of such pertur- 
bations at a given a is achieved by the perturbation 
which is reaching its virial size at a, and is equal to 
8 v (a), given in line 3 of Table II. The applicabil- 
ity domain of the second branch of equation C15 is 
then fe(a) < 5 < fe(a) and the conversion relation 
in this case is shown in the 3rd line of table III. 

• Perturbations which have reached their virial size 
but have not yet reached their designated collapse 



time, acoii(A0> need to be treated separately, since 
the spherical collapse model fails (does not agree 
with the physical picture we would like to describe, 
although it is still formally applicable) for radii 
smaller than the virial radius. Since a realistic de- 
scription of the microphysical dissipation processes 
which lead to virialization is far beyond the scope 
of this analytical calculation, we will adopt a pre- 
scription which is driven by mathematical simplic- 
ity. We will assume that for S v (a) < 5 < 8 c (a) the 
conversion relation So(a,5) has the simplest poly- 
nomial form which satisfies the following physically 
motivated boundary conditions: 

— The extrapolated overdensity is continuous 
and smooth at S v , so <5 (a,<5 v ) = 8o, v and 

d&o/dS = 8' Qv as given by the appropri- 

<5 V 

ate relations of the previous branch (3rd line 
of tables III and IV correspondingly). 

— After time a co ii(/x) the radius of the perturba- 
tion remains constant and equal to the virial 
radius so changes in the (true) overdensity are 
only due to the increase of the scale factor of 
the background universe. This then implies 
that So(a 1 S c ) — <5q,c given by equation C29 and 



OSo/dS 



(d~8 /d^ c )(d»/d8\ Sc ) = ~8' c 



O.c 



given in Table III line 4. 



The conversion relation in this case is shown in 
Table III line 4, while its applicability domain is 
S v (a) < S < S c (a), with S c (a) given in Table II 
line 4. Past their collapse time, perturbations are 
treated as virialized objects without substructure 
and are not relevant as "local environment" of other 
objects for the purposes of our double distribution 
calculation, hence it is not necessary to have a con- 
version relation of 5 > 8 c (a). The calculation would 
be simplified (this last branch would be unneces- 
sary) if we chose to regard perturbations as virial- 
ized objects after the moment they reached their 
virial size after turnaround, at time a v . However, 
we will retain the usual assumption that objects 
virilize at time a co ii for consistency with existing 
Press-Schechter calculations. 



3. Critical extrapolated overdensity for collapse, 

<5o,c(a) 

We need to find the critical 8o tC (a) for collapse if the 
field is linearly extrapolated to the present epoch, i.e. the 
value the linearly extrapolated to the present overdensity 
must have, for a structure to have collapsed at universe 
scale factor a. This, from equation C28, will be 



3A[(2u)W] 1 + Mc (q) 



(C29) 
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TABLE II: Applicability limits for different branches of the spherical evolution solution 



Limit 



Expression 



Auxiliary Relations 



Density contrast of 
Eddington overdensity, 5Ed(o) 



5 Ed (a) 



1 sinh _1 Vcuo 3 - Vi (r Ed , 2) = 



Density contrast of overdensity 
turning around at a, 5t*(a) 



5 ta (a) = Lua 3 n ta ,(a) - 1 



sinh 1 Voja 3 — Vi(l,p ta ) = 



Density contrast of overdensity reaching 
its virial size at a, 5 v (a) 



S v (a) 



«P,v[fv(n)] 



sinh 1 s/oja? — 2Vi(l,/U v ) + Vi [r v (/i v ),pi v ] = 
r v (/j) = ap,v(/x)/a P ,ta(p) 



Density contrast of virialized overdensity 
formally collapsing to a point at a, <5 c (a) 



5 c {a) = 



ip,v[Mo(a)] 



- 1 



sinh" 1 \/wa3 - 2Vi (1, Mc) = 



TABLE III: Different branches of conversion relation 5o(a, 5) for an f2 m + Qa = 1 universe. 



Branch 



8o(5, a ) 



Auxiliary relations 



— 1 < S < 8 Ed So(a, 5) = 



3A[(2 W )V3] 


1 — zu(a,5) 


2 l/3 


Ma,<S)] 2 /3 


3A[(2^) 1 / 3 ] 


l+/j,(a,6) 


2 l/3 


\„(a,S)] 2 ' 3 


3A[(2^) 1 / 3 ] 


l+/»(o,«) 


2 l/3 


[ M (a,«)]2/3 



sinh 1 s/oja 3 — Hi 



( X 1 / 3 

a {—) > w 



= 



fed < <5 < (5 t a 



<S (a,<5) = 



sinh 1 



Vi 



1/3 

T+J / ' A 4 



= 



<5ta < 5 < <5v 



5 (a,5) 



sinh 



/W -2Vi(1,m) + Vi 



, 1/3 
1+5 / 'A* 



<5 V < <5 < 5 C 



<5o(a,5) = 5o,v + 5q >v (5 — 5 V ) 

+ — — ; — («c-«v)^ : — ; : — ^ ~ K > 

I ( l5 0,c + l5 0,v)( ,S t:-4v)-2( 1 5 0iC -5o, v ) , x x ,3 

+ (-Sc-M* (O-OvJ 



5o,v = So(a, 5 V ), given in this Table line 3 
5o,v = I , given in Table IV line 3 

1 <5 V 

So, c = So(a,8 c ), given by equation C29 



1/31 



, . „„,2/3„2 „2/3_ 
^ = W U C ~ (2Ml) 1/3 (ic + D 

/ic(ffl) given in Table II line 4. 



l-Mc) 



where again /i c (ct) is given by Table II line4. 

The dependence of $o, c on a can also be expressed in 
terms of the linear growth factor, D(a), as was the case 
for the O m = 1 universe. The conversion relation between 
<5(a) and 5 a (the linear-theory result for the density con- 
trast at time a) is independent of a. In other words, as 
long as 5 and 5 both refer to the same time, knowledge 
of the one uniquely defines the other, independently of 
the actual time at which they are both evaluated. As 
5(a) — > 00, 5 a — > 5 C , the linear-theory density contrast at 
the time of collapse (given by eq. C29 for a = 1). There- 
fore 5 C is the same for perturbations of all curvatures, 
and, using equation C26, we can write 



<W(a C on) 



^ D(ao) 

° D(Ocoll) 



(C30) 



4. dS /d8\ a 

In addition to the relation between <5 and <5o, we will 
also need the derivative d8o/d5\ a in order to convert be- 
tween true and extrapolated overdensity differentials in 
equation (16). The calculation is similar as in the case 
of an fl m = 1 universe, and the results are summarized 
in table IV. 

In table IV, H.2(r,m) is the hyperbolic vacuum integral 
of the second kind, and V2 (r, ji) is the incomplete vacuum 



integral of the second kind, defined in appendix D. 



APPENDIX D: VACUUM INTEGRALS 

Due to the central importance and frequent use of vac- 
uum integrals in this calculation, we have generated nu- 
merical functions based on combinations of tabulated val- 
ues and asymptotic approximations which evaluate each 
of the vacuum integrals in a small fraction of the time 
that would be required for quadrature, and with an accu- 
racy better than 0.5% throughout their domains. These 
functions are publicly available at http : //www . astro . 
uiuc . edu/~bdf ields/DD. 



1. The incomplete vacuum integral of the first 
kind Vi 

a. Definition 

We define the incomplete vacuum integral of the first 
kind as 



Vi(r. 



- X)(-X 2 - X + Jbj 

with domain < r < 1 and /i > 2. 



(Dl) 
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TABLE IV: Different branches of derivative d8o/dS\ a for an fi m + ft a = 1 universe. 



Branch 



Auxiliary Function 



d8o/d8 (a, <5) = Definitions of Additional Functions 



-K8<S Ed $! =sinh- 1 ^^-'Hl(r, ro ) (ffi) ^ 



a 3/2 



2(1+5)3/2^/^ 

/ \l/3 
/ row \ __ 

a {T+l) 



n_i 

as ' sa^uj) 1 / 3 ] n2 



sinh 1 \fuia 3 — Hi 



fed < S < fe $2 = sinh 1 Vua 3 - Vi (r, u) 



2/3(1 + 5)1/3 ■ 
1/3 



a<s 



„3/2 



2 ( 1 + ^ /2 ,/F&)- u 2/3 ( ( ^ ) ) 1/ 3 - 



Si 3A[(2t,;) 1 /3 

sinh -1 V<^a 3 — Vi 



V 2 



1+5 



1/3 



\ 1/3 



fe < <5 < fe 



$3 = sinh 1 yu. 



-2Vi(l,Ai)+V(r,Ai) 



(> r- 

95 



V 8S J 



9*3 



a 3 / 2 





2(l + 5)3/2 v 


„ 3 




a(/.+ l) 








fl 2/3 (1+5) l/3 i 




d<S>3 


(2^)1/3 


I- 


a 




1/3 




8S 


3^1(2^)1/3] ) 







6y dVi(l,M) 



sinh" 1 Vua? _ 2V1 (1, tt) + Vi 



/ \ 1/3 



5 °.v + 2 (5^K? ' ( S - Sv > With So,c, 5o,v 5(',,c, *0,v, 

, (^,c+^,v)(^-^v)-2(io,c-^o,v) ^ x , 2 as in Table III line 4 
(i c -f v )3 [d-o v ) 



fe < 5 < S c 



+3- 



b. Properties 

Physically, Vi (r, /i) is proportional to the time re- 
quired by a perturbation of normalized curvature pa- 
rameter k/uj 1 / 3 = (ti + l)/ii 2 / 3 to achieve a size a p = 
rap^zin/u 1 / 3 ) before turnaround. Its asymptotic behav- 
ior for r < 1 is 



Vi(r, M ) 



r«l 1 

\/m' 



,.3/2 



(D2) 



while for \i ^> 1 it is 



Vi(r, M ) 



(U»l 1 



— \/r{l — r) — sin 1 \f\ 



(D3) 

In the case of the Eddington perturbation (fi = 2), we 
can derive a closed- form expression for Vi : 



Vi(r,2) 



y/xdx 



2 J (l-x)Vx + 2 
V3 
2 



l n i±^_2V3sinh- 1 



+ ln 



W 

2V3 + v/3r + 3V2T7 
2V3 - v/3r + 3V2T7 



(D4) 



When r = 1, the value of Vi(l, /i) is the complete vacuum 
integral of the first kind, which is a function of /x alone. 
Physically the complete vacuum integral of the first kind 
is proportional to the time required for a perturbation of 



curvature parametrized by \x to reach turnaround. The 
derivative of Vi(l,/z) appears in the calculation of the 
derivative dd /dd, in the 3rd line of Table IV, and it is 



— v a ) = -- C ^ dx 



(D5) 



2. The hyperbolic vacuum integral of the first kind 

Hi 

a. Definition 

We define the hyperbolic vacuum integral of the first 
kind as 



*■<'■"> -5jf:7iP 



y/xdx 



(D6) 



/0 + x){x 2 - X + w) ' 

with domain < r < oo and w > 1/4. 

b. Properties 

Physically, Hi(r, w) = is proportional to the time re- 
quired by a perturbation of normalized curvature pa- 
rameter k/oW 3 = (1 — vj)/w 2 ^ 3 to achieve a size a p = 
?"%>,r(k/w 1/ ' 3 ). Its asymptotic behavior for r <C 1 is 



(D7) 
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while for r > 1 it is 

Hi (r, w) r i 1 C{w) + | In (Wr 2 - r + w + 2r - l) 

(D8) 

where C{w) is a function dependent only on w. In the 
case of a flat (w = 1) perturbation, Hi can be integrated 
immediately to give 

Wi(r,l) = § r^^==sinh- 1 V^. (D9) 
2 Jo Vr + 1 

3. The incomplete vacuum integral of the second 
kind V 2 

a. Definition 



4. The hyperbolic vacuum integral of the second 
kind Hi 

a. Definition 

We define the hyperbolic vacuum integral of the second 
kind as 

*/ / \ 3 r x 3 / 2 dx 

H2{r ' W) = -J {l + X )W(a?-x + «,)*/» (DU) 

and its domain is that of H\(r, w). 

b. Properties 



We define the incomplete vacuum integral of the second The hyperbolic vacuum integral of the second kind is 
kind as related to Hi(r, w) through 



3 [ r x 3 / 2 dx 

V2(r ' M) = 4/ (i-xr/H-x 2 -x+^ > (D10) o 



with domain same as for Vi (r, . 



d{n/^/ 3 ) 



Hi(r,w) = vo 2 l 3 H 2 (r,vo), (D15) 



b. Properties 



with 



The incomplete vacuum integral of the second kind is 
related to Vi (r, /tz) through 



d 



din/to 1 / 3 ) 



— Vl (r^) = ^ 3 V 2 (r^) , (Dll) 



with 



k/^ 3 = (1 + ^/^ 3 . 



(D12) 



In the case of the Eddington perturbation (fi = 2), we 
can derive closed- form expressions for V 2 : 

, . 3 r x 3 ' 2 dx 

V2(r ' 2) = 4j (l-x) 3 (x + 2) 3 / 2 



V3 
72 

+ log 



/37(2-3r + 4r 2 ) 
V2T7(1 -rf~ 

1-r 



1 + 2r + v/3r(2 + r) 



(D13) 



K/u 1 ' 3 = (l + m)/w 2 ' 3 . 



(D16) 



In the case of a flat (tn = 1) perturbation, takes the 
form 



H 2 = 



2 4/3 



f ( 



u 3 + 2 



3/2 



(D17) 



which is the integral entering the linear growth factor in 
the Q m + Cl\ = 1 universe. Hence, the linear growth 
factor function A(x) can be written as 



A[x) = 24/3 ^ + 2 ) 1/2 H 2 (2^/3 X; 1} . (D18) 



